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1. Introduction 

Landauer-Biittiker (LB) scattering theory pQ marks a powerful approach to multi- 
terminal transport in conventional (Fermi liquid) mesoscopic devices, but is well-known 
to be inapplicable to strongly correlated electron systems. Progress along this direction 
has recently been reported for interacting metallic ID systems, where a good description 
in terms of a Luttinger liquid (LL) |2 a is often available. Single-walled carbon nanotubes 
(SWNTs) represent a prime example for such a system [HI HIE]- F° r a two-terminal setup, 
the generalization of the Landauer formalism to a LL (containing an impurity) coupled 
to noninteracting leads has been formulated and solved previously [HI Ej . For related 
work, see Refs. |Hl El HOI HH H21 HS1 H31 [1^1- 

Recently, these ideas were generalized to also treat iV-terminal starlike 
LL/nanotube networks (N > 2) 16 . Such a transport theory is then able to make 
concrete predictions for experiments, e.g. on carbon nanotube Y junctions. First steps to 
the production and experimental analysis of SWNT networks are under way. Template- 
based chemical vapor deposition [T7] and electron beam welding methods jTHj allow 
to fabricate and contact multi-terminal nanotube junctions. Recently, an intrinsically 
nonlinear current- voltage characteristics of multi-walled nanotube Y junctions has been 
reported ^H]- Furthermore, four-terminal devices have been realized by several groups 
using two crossed NTs [201 EI]- Such nanotube-based (or related) setups are of growing 
interest to many theorists working in the field of interacting mesoscopic systems. To 
mention some recent work, in Ref. [22] a detailed renormalization group analysis of a 
similar model as ours was performed. Other authors used perturbation theory in the 
hopping [23121] or in the interaction strength [22]. Furthermore, the conductance of a 
junction of three Luttinger liquids enclosing a magnetic flux was investigated [2E] . The 
cited activity underlines the need for a generalization of the LB approach to correlated 
systems. With this work we hope to contribute a significant step towards this aim. 
Some of our results have been published before ^B] , but most of the material is new. 

The paper is organized as follows. In Sec.[2]we describe our model for the interacting 
multi-terminal nanotube device, and outline its general solution in the symmetric case. 
As the step from N = 3 to iV > 3 is straightforward JH], we concentrate on nanotube Y 
junctions (N = 3) in the following sections. In Sec. E] the nonlinear conductance matrix 
of a symmetric Y junction is discussed. The case of asymmetric junctions is then studied 
in Sec. HI Finally, we conclude and point out some open questions in Sec. El 

2. General model 

We shall first discuss the simplest case of a spinless single- channel Luttinger liquid. In 
fact, in the first part fSec. 12.1]) . we even ignore the interactions and briefly review some 
aspects of LB theory relevant to our purposes. The case of a symmetric junction is then 
treated by various boundary conditions for the interacting system in Sec. 12 .'2\ and the 
extensions necessary to treat nanotubes are discussed in Sec. 12.31 These include the 




Figure 1. Basic setup: N nanotubes form a starlike junction (shown for N = 3). 
Each nanotube is supposed to be adiabatically contacted, with chemical potential fii 
of the respective reservoir. The cross indicates an additional impurity. 



electronic spin and an additional flavor degree of freedom due to the presence of two K 
points in the first Brillouin zone jlj. The basic setup studied here is shown schematically 
in Fig. [TJ 

2.1. Fermi liquid case 

When one ignores the Coulomb interaction among electrons, transport in a starlike 
setup of N nanotubes can be described by the scattering matrix approach, employing a 
unitary S matrix with entries Sy determining the transmission coefficients = \sij\ 2 
of the system. The current through nanotube i may we written as 

h = T JZ dE ( [1 ~ Tii ^ E ^ E - - £ T v( E )f( E - /*)}. 

where f(E) is the Fermi function (we consider only a single mode here, while a true 
SWNT has four). Unitarity of the S matrix implies (1) the Kirchhoff rule (charge 
conservation), X^^V/ = 1, and (2) gauge invariance under a uniform potential shift in 
all the reservoirs, J2j Tij = 1. When all transmission coefficients are energy- independent, 
the nonlinear conductance matrix normalized to e 2 /h reads 

_ hdh_ _ m 

Without applied magnetic field, time-reversal symmetry implies = Tji. For concrete 
calculations of the S matrix for nanotube Y junctions, we refer to Refs. (23 EU 
The N x N matrix S relates the outgoing and incoming scattering states at the node, 

* £ (0) = S* R {0). (2) 

To describe the N (for now semi-infinite) tubes, we use ^(x) = (ipx, . . . , ip^) — + 
with x < 0, which contains the outgoing (left-moving, L) components as well as the 
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incoming (right-moving, R) states, with the junction at x = 0. Unfortunately, the 
boundary condition (J2J) is exceedingly difficult to handle for correlated electrons. 

To make progress, we shall first consider only symmetric junctions. A wide class of 
such junctions can be parameterized by a symmetric S matrix of the form 



' N+iX x N+iX N+iX > 

2 2 1 2 



s 



N+iX N+iX " N+iX 



(3) 



V _2_ _2_ . . . _2 1 / 

where A > serves as arbitrary parameter. This is the most general S matrix that 
implies wavefunction matching at the node [3U| . 

Vi(0) = ... = V*(0), (4) 
for the components of \I/ = ^> l + The S matrix (jSJ) leads to 

T. - ( iV - 2 ) 2 + A2 T . - 4 ^ 

iV 2 + A 2 ' * J ' iV 2 + A 2 ' 1 } 

which, for A = 0, describes a maximally transmitting symmetric junction, T# = 
(N — 2) 2 /N 2 and T^j = 4/N 2 . On the other hand, for A — > oo, a very bad junction 
forms, Tjj = 1 and Tj^j = 0. For arbitrary A, the conductance matrix then follows from 
Eq. ((TJ). It is worth noting that no symmetric starlike junction can act as an ideal beam 
splitter, since backscattering at the node is unavoidable. The only possibility to have 
an ideal beam splitter, T# = 0, arises in asymmetric junctions 

Since phases in ^(0) can be gauged away, a density matching condition arises for 
any symmetric junction described by Eq. (J3J), 

Pl (0) = ...= pjv (0), (6) 

where Pi(0) denotes the total electronic density in tube i at the node. The condition (pj) 
is crucial in our approach for interacting systems, as it avoids the explicit wavefunction 
matching. The density in Eq. © refers to the sum of the slow (long wavelength) part, 

= *JM*fl,< + *1<*L,<> and the Friedel oscillation, 5 Pl (x) = ^ R ^ LA + ^I^^i- 
Since the parameter A will affect the Friedel oscillation, there is a hidden dependence 
on A not directly visible in Eq. (jOJ. 

2.2. Luttinger liquid case 

Next we include the bulk interaction in the tubes, but still restrict ourselves to a 
symmetric junction. We assume that all nanotubes can be characterized by the same 
interaction strength, and for now still consider just the spinless single- channel version. 
Following Ref. we may neglect backscattering and therefore obtain a Luttinger liquid, 
where interactions are encoded in the standard Luttinger parameter g < 1, with g = 1 
referring to the noninteracting limit. Under bosonization j2], the Hamiltonian is (we 
put h = 1) 



2 hi-** 



tf + - 2 {d x 6if 



(7) 
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where Hi(x) is the canonical momentum to the standard boson field 6i(x) for tube i, 
and vp is the Fermi velocity. Similar to the two-terminal case, adiabatically coupled 
external voltage reservoirs (at x ~ — L) lead to radiative boundary conditions for the 
long- wavelength part of the chiral electron densities [B], 

(g- 2 + l)pl R {-L) + GT 2 - l)p i>L (-L) = m/irvp. (8) 

These conditions only depend on the injected currents, which are independent of any 
backscattering happening at the node or within the nanotube. To fulfill Eq. (JH|) . the 
chiral densities must combine to give 

h = «>f(pIr ~ pIl) = ^[vi~X f « U ^ ' ^ 

P° = PU + PIl = ^U + Y, > ( 1Q ) 

which defines effective transmission matrix elements Tij(g, {pk}, ksT). The matrix T 
should not be confused with the LB transmission matrix T of the noninteracting system. 
Albeit T — T for g — 1, the effective transmission matrix T cannot be interpreted as 
single-electron transmission matrix for g < 1. The applied voltages in Eqs. Q and (fTU|) 
are defined as 

eUi = fii-p,, fi = J2pi/ N - (11) 

i 

Under this definition, gauge invariance is automatically fulfilled if the transport 
properties of the system only depend on the U{ and not explicitly on the {fik}- Charge 
conservation (Kirchhoff 's rule) can be incorporated by imposing 

N 

£^ = 1- (12) 

1=1 

Finally, the density matching conditions (JHJ) have to be obeyed. 

Due to the pinning amplitude 5pi(0), i.e. the value of the Friedel oscillation at 
the node, Pi(0) will deviate from p° in Eq. (|ID)). To compute the pinning amplitude, 
we first note that the Friedel oscillation 5pi(x) in tube % arises due to interference 
of the incoming right-movers and the left-movers that are backscattered at the node. 
Left-movers transmitted from the other iV — 1 tubes into tube i cannot interfere with 
the incoming right-movers and will therefore not contribute to 5pi(x). To verify this 
argument explicitly, it suffices to analyze the corresponding two-terminal case with 
one impurity at x — 0. Without loss of generality, let us take /c-independent bare 
transmission amplitude t for the impurity. Taking k > 0, there is a right-scatterer 
eigenstate 4>k(x), 

if e lhx + re~ thx (x < 0) 



\ te ikx (x > 0) ' 

where \t\ 2 + |r| 2 = 1. Similarly, the left-scatterer is 0_fc(x) = (/>k(—x). We now show 
that there is no contribution to Sp(x) from interfering right- and left-scatterers. While 
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this statement is immediately seen for noninteracting electrons, it can also be confirmed 
for g < 1. To that purpose, let us expand the electron operator in terms of the right- 
and left-scatterers instead of the standard R/L movers, 

^( x ) ~ a k<t>k{x) + b k <f>-. k (x). (13) 

fc>0 

Of course, the kinetic part of the Hamiltonian decouples in the Fermi operators at 
and bk- Moreover, when inserting Eq. (|13|) into the interaction Hamiltonian, no terms 
containing only one or bk can appear because the right- and left-scattering states 
4>k{x) and are orthonormal. As a consequence, in the expectation value of the 

density operator, there can be no terms involving (a\b k ) or (p\a k ). Hence there is no 
interference between scattering states originating in different reservoirs, a finding that 
directly translates to the TV-terminal setup under study here. 

As a consequence of the above discussion, 5pi(0) must be identical to the pinning 
amplitude in a related two-terminal setup defined by an impurity whose bare reflection 
probability is determined by Using this mapping, we can thus compute the pinning 
amplitude from 

(«Jpi(0)> = -gfeVi/nvp, (14) 

where the four-terminal voltage V% = Vi(Ui) has to be determined through a self- 
consistency relation, see below. The validity of Eq. ()14j) can be seen as follows. The 
long-wavelength density fluctuation in the equivalent two-terminal setup, due to an 
applied voltage Ui as given by Eq. (JTTJ), is [E] 

(p°(x)) = -g 2 eVi sgn(x)/2nv F . 

With respect to a density corresponding to the chemical potential £L, taking x < 0, we 
thus have p° = g 2 eVi(Ui)/irvF- For the two-terminal problem, antisymmetry of (p{x)) 
under U — > —U implies that the net change in density due to U is zero at x = 0~. 
Therefore p° must be cancelled by the pinning amplitude, see Eq. ((HJ). Note that this 
reasoning does not assume the validity of the density matching conditions ©• 

For g = 1/2, the self-consistency relation determining V(U) can be obtained using 
refermionization, and takes a rather simple form. Results for this special interaction 
strength have been discussed in Ref. ^Hl- Using integrability concepts, however, the 
respective self-consistency relation can also be derived and solved for arbitrary g [Zj. 
Let us write down the necessary equations for g = 1/p (with integer p) , keeping in mind 
that arbitrary values of g can also be treated. The bare reflection probability of the 
impurity determines the energy scale k B T B oc exp9 B , see below, and the elementary 
excitations scattered one-by-one at the impurity are kinks, antikinks, and p — 2 neutral 
breathers. The two-terminal current I then follows as an integral over rapidities 9, 

I = J \T ++ \ 2 (a + -a_)d9, (15) 

where |T ++ | 2 = 1/{1 + exp[— 2(p — 1)(9 — 6b)]}- The occupation probabilities are 
o~± = nf± for the kinks/antikinks, where 2irn = de/d6 with the eigenenergy e of 
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the kink/antikink. Explicit expressions for the eigenenergies can be found in Ref. [7j. 
Furthermore, the filling fractions read f±(9) = 1/(1 + exp[(e(6*) =F W/2)/T]), where 
±W/2 is an effective chemical potential for the kink/antikink. The parameter W has 
to be determined self-consistently through 

J (|T ++ | 2 + p|T+_| 2 ) {a+ - o_)d6 = geU/2nv F , (16) 

where U is the applied voltage and |T + _| 2 = 1 — |T ++ | 2 . Since (e 2 /2tt)(U — W) = (1— p)I, 
see Ref. 0, the solution of Eqs. JTSJ and then yields W = W(g, k B T, U), and finally 
the four-terminal voltage follows from V = W — 2irl/ge 2 . 

In these equations, Tb appears as effective impurity strength of the related two- 
terminal problem. With Ta given in Eq. (jSJ), an effective dimensionless impurity strength 
vq can be determined by examining the noninteracting limit, see below. Given vq, the 
scale Tb then follows as 

k B T B /D = c 9 vl /{1 - 9 \ (17) 

where D is the electronic bandwidth and c g a numerical prefactor of order unity whose 
precise value is given in Ref. [Zj. Using this correspondence, the complete transport 
problem for a symmetric junction described by the bare S matrix (jSJ) can be solved 
explicitely. Using Eq. (JHJ), the density matching conditions © read (l + Tjj)Uj — 2Vj = 
TkjUj for all pairs k ^ j. Exploiting Kirchhoff's rule JEJ), this set of equations admits 
only the solution 

~ _2-N 2(N-m ~ _l_2K 

u ~ N + NUi ' Mi ~ N NU; 1 ' 

Evidently, the depend only on the applied voltages Ui but not separately on 
the chemical potentials /ij. Thus gauge invariance is automatically fulfilled. The 
conductance matrix then follows immediately from Eq. ©. In Sec. El we will discuss 
this solution in some detail for the case N = 3 (nanotube Y junction). It then only 
remains to fix the impurity scale Vq employed in determining Vi (£/*). 

As the impurity strength parameter vq can be obtained in the noninteracting limit, 
a convenient way to fix vq is as follows. For g — 1, as pointed out above, Ta should 
simply reproduce the LB transmission coefficient Ta specified in Eq. (J5J). In addition, in 
the noninteracting limit, the ratio VI/ Ui is just the bare reflection coefficient R of the 
impurity in the two-terminal problem [6 . The latter can be expressed in terms of vq 
according to 

R- «» 



[i+vinr 

where Vq = 2 corresponds to the unitary limit; larger values for vq lead to unphysical 
results. Combining above equations then leads to 

u (A, N) = 2{VN 2 + A 2 - V2N)/^jN{N - 2) + A 2 . (19) 

For a very bad junction, A — > oo, we are indeed in the unitary limit, Vq = 2. For a 
perfect junction with A = 0, however, v = 2(y/N — \/2)/y/N — 2. In particular, for 
iV = 3, this value is v ~ 0.63567. 
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2.3. Extensions 

Real SWNTs are characterized by four channels, resulting from the electron spin and the 
additional flavor degree of freedom. Including these degrees of freedoms as additional 
indices (a, a), Eq. (J2J) now involves a AN x 4N dimensional S matrix, and the electron 
state \& = tyft + has 4N entries ?/w(aO- Assuming that at the node neither 
flavor mixing nor spin- flip processes are present, the previous analysis can largely be 
carried over. Since the S matrix factorizes, the Kirchhoff rule as well as the density 
matching conditions hold separately for each component (era). The Hamiltonian (|7j) 
then corresponds to a four-channel boson model 0] in terms of symmetric/antisymmetric 
charge/spin fields. The interaction parameter in the symmetric charge mode is then 
g c+ < 1, while the neutral modes are essentially at the noninteracting value. The 
radiative boundary conditions (j8|) are only imposed in the symmetric charge sector, 
and our above discussion then directly leads to a very similar expression as in Eq. Q. 
However, when computing the effective transmission matrix T, the four-terminal voltage 
parameter Vi is now nontrivially affected by the additional degrees of freedom. This 
is due to the fact that backscattering at the node couples the four boson modes 
diagonalizing the Hamiltonian. At the moment, we are not aware of an exact solution 
to this problem. Nevertheless, two remarks can be made even in this situation. First, 
the prefactor in the current © is changed from e 2 /h to 4e 2 //i, reflecting the presence of 
four channels. Second, we expect that the typical power law exponents to lowest order 
are modified according to g — > (g c+ + 3)/4. Therefore, the qualitative physics seems to 
be captured already by the one-channel version discussed further on. However, we wish 
to stress that the finite value for vq obtained under the mapping does not seem to allow 
for simple perturbative access to the full four-channel problem. 

3. Conductance of a symmetric Y junction 

In this section we discuss the conductance defined by Gy = (h/e)dli/dfij for the case 
of a symmetric Y junction (N = 3), based on the solution of the general problem in 
Sec. 12.21 Using Eq. (j!8)) gives the nonlinear conductance matrix in the form 



where k ^ i ^ j in the second equation. Eq. (|2()jl holds for any interaction strength g. 
To obtain explicit results, one only needs to compute V(U) for the corresponding two- 
terminal problem. We shall show explicit results now for several values of A, the special 
interaction parameter g = 1/3, and applied voltages \i\ = E F + elf, ^2 = ^3 = E F . 
In Fig. |2J the linear conductance is shown, and in Fig. El the nonlinear conductance is 
plotted for k B T/D = 0.01. 




(20) 



le-K 



1.01 



0.1 



k B T/D 



Figure 2. Temperature dependence of the linear conductance coefficient G\2 for 
g = 1/3 on a double-logarithmic scale, for several different values of the parameter A 
□ 




Figure 3. Voltage dependence of the nonlinear conductance Gn (double-logarithmic 
scale) at ksT/D = 0.01, for several values of A. 
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As illustrated in Figs. 121 and El for G 12 at g — 1/3, at low energies the system always 
flows to a fixed point corresponding to disconnected nanotubes, 



To the best of our knowledge, all theoretical work on this problem [To^ 12*2*1 12*3*1 l2~H 12*5^ l2*B] 
agrees that for symmetric junctions, Eq. (}2*Tj) represents the only stable fixed point in 
the case of repulsive electron-electron interaction. This fixed point represents an isolated 
node weakly coupled to N broken-up quantum wires. It is then of interest to study the 
energy dependence of the Gij at low energies, where a typical LL (^-dependent power law 
behavior arises, see also Ref. [23] . The leading corrections to Eq. (}2*Tj) can therefore be 
formulated in terms of the well-known end tunneling density of states, with associated 
power law exponents given in textbooks 

4. Conductance of an asymmetric Y junction 

Next, we construct more general S matrices starting from the solution for the symmetric 
situation discussed in the two previous sections, focussing for clarity on the case N = 3. 
The basic idea to make the system asymmetric is to add impurities of dimensionless 
strength Vi to tube i. To avoid the appearance of spurious resonances, these are taken 
sufficiently close to the node, at, say, x ~ —1/kp- This modelling of an asymmetric 
junction still allows for g < 1, e.g. by using perturbation theory in the around the 
above solution for a symmetric junction (vi = 0). The prototypical situation for exactly 
one impurity, v i > but v 2 = v 3 = 0, is illustrated in Fig. [TJ Since the impurities are 
displaced away from x = 0, the junction itself is characterized by the symmetric S matrix 
(JHJ), and therefore the density matching conditions (jBJ) can still be used. Nevertheless, 
effectively this procedure allows to generate asymmetric S matrices. 

The free real-time action then reads in the framework of bosonization [2j 



where the bosonic phase field 6i(x, t) is decomposed according to 9i = 9\ + 9f . Here, Q\ 
obeys homogeneous {Ui = 0) boundary conditions (jH)), and 9f is a particular solution to 
the equations of motion obeying Eq. (jSJ). Written in terms of the boson fields, one has 
I®(x,t) = (e/y/n)d t (0i(x, t)) and p°(x,t) = (1/ y/7r)d x (9i(x, t)), and one checks easily 
that the equations of motion as well as Eq. (jHJ) are satisfied by 



where If and p° are given by Eqs. © and (|TU|). respectively. Here we assume that 
applied voltages are chosen such that fl = 0. 

Adding an impurity of strength A x = Dvi/tc (with bandwidth D) into tube 1 then 
yields the action contribution 



(21) 



s [0i] f dx i dt \ (<W 2 - -2 (W 2 

£* ■ -i J —L J Q 



(22) 




(23) 




(24) 



Nanotube Y junctions 



11 



^From now on, L — > oo to allow for analytical progress. Inserting Eq. (|23|) into Eq. (J22|) . 
we obtain 

s Q [o t ] = s m - vW£(Mo)/<? 2 ) J dt e}(o,t). 



i=i 



The second term can be dropped because of Eq. (JUJ) and Kirchhoff 's rule, which implies 
that J2i 9^(0, t) equals an arbitrary constant (which can be put to zero). The current Jj 
is then given by 



.-,2 



(25) 



4-1. Weak asymmetry: Perturbation theory 



The impurity strength Ai measures the degree of asymmetry of the junction. First, 
we treat the case of weak asymmetry under a perturbative calculation in Ai, using the 
Keldysh technique. In this subsection, we write 9 h —>■ 9 and restrict ourselves to T = 0. 
Then the current ()25|) follows from 

(dA(o,t)) = \jl\ V9 i e iS ^5 (l>(M)) dA(0,i), (26) 

i=l \ i / 

where S[9i\ = So[9i\ + Si. The 5-function reflecting Kirchhoff's rule is enforced 
using a Lagrange multiplier field j(t), yielding the additional action contribution 
S 1 = J dt r f(t)J2i9i(0,t). Switching to Fourier space, the action S[9i,j] = So[9i] + 5 7 
reads 

£ % [i g - ?) m + . 

When switching to fc-space, a principal value contribution has been dropped, which is 
justified in the low-energy regime of interest here • We then introduce rotated phase 
fields, Q± = 9 2 ± 9 3 , with conjugate momenta U± = (n 2 ± n 3 )/2, and subsequently 
integrate out the fields 9\ and 7 by standard Gaussian integration. This procedure 
leads to the (impurity- free) effective action So^+j + SId [<?-]> where 

*w-f/s5(5-S)^- (27) 

with v + = 4:Vp/3,g + = 4g/3,v_ = 4vp, and g_ = 4g. Fourier transformation back to 
real space gives 



S [9 ± ] = ^ r dxdt 
2 J —00 



\(d t 9 ± f - \{d x 9 ± f 
v± g± 



(28) 



where the rotated boson fields 9±(x) live on the full line. As we shall see below, only a 
certain linear combination of the right- and left-moving components of 9 + couples to the 
impurity, and hence the extension to the positive half-line (x > 0) does not introduce 
spurious effects while simplifying the subsequent analysis considerably. 
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Expectation values of the form (|26|) are then calculated by first evaluating 

(d t 9±{t)) = \ [ V9 + V9_e i{s[e+]+So[e - ]) d t e ± (t), (29) 
where the impurity in tube 1 acts only in the + sector, 

S[9 + ] = S [9 + ] + Ai / dtcos[-V^r9+(0,t) + 2irl%t/e]. (30) 

Since (d t 9-(t)) = 0, we are left with the relations (d t 9i) = —(d t 9 + ) and (d t 9 2 ,3) = 
(d t 9 + )/2. Equation (|2~9")l can now be evaluated perturbatively in Ai, and Eq. (|23j) then 
gives [HI] 

I 1 = I?-6I Al , I 2>3 = I% )3 + 6I M /2, (31) 

where If denotes the currents computed in Sec. 12.21 in the absence of the impurity, and 
the lowest-order perturbative correction due to the impurity is 

6I Al = eg + (k\/D) sin(27r^ + )r(l - 2^ + )sgn(/ 1 )(27r|/ 1 °|/e J D) 2 ^- 1 . (32) 

To treat more general S matrices, one can add a second (weak) impurity A2, say, 
into tube 2. Then, within order Af, again correction ()32j) arises, the correction of order 
AiA 2 vanishes, and corrections of order A 2 , may be calculated in the same way as above. 
This gives 

= -e( 7± (A2/2 J D)sin(2vr( 7+ )r(l-2( 7+ )sgn(/ 2 )(27r|/ 2 |/e J D) 29+ - 1 , (33) 

where we replaced most of the g_ that actually appear in Eq. (|3*3*|) by g + , utilizing 

g + + g~ = 4g + , and the currents Jj are 

h, 3 = I° 2i3 + (6I Al + 8I£ 2 ±5I^)/2. 

Evidently, a wide class of asymmetric Y junctions is accessible by this construction 
method even for interacting SWNTs. 

4-2. Arbitrary asymmetry: Refermionization 

In the case of just one impurity (Ax), it is possible to solve the full transport problem 
in a nonperturbative manner at the special interaction strength g = 3/8. This point 
corresponds to the Toulouse limit g + = 1/2 in the rotated picture used in Sec. 14. 1| 
and allows to use refermionization e 2j. In fact, one can solve this transport problem for 
arbitrary g using integrability methods, see Ref. [Zj, but for clarity, in this section we 
only discuss the simplest non-trivial case, g = 3/8, which admits the use of standard 
techniques. We only have to study the 9 + sector, as the 6L field decouples from the 
impurity and can be ignored. 

To sketch the g = 3/8 solution, let us first define new chiral boson fields (<fi+(x) is 
the conjugate field, 11+ = d x (p + ) 

Vr/l{x) = \/Tr(y/g^<f> + (x) ± 9+(x)/ ^ ), 
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which obey the algebra [(p p (x),ip p i(x')] = —iir5 pp isgn(x — x'). Next we switch to linear 
combinations of these chiral fields, (p± = (Lp R =p ip L )/\/2, with associated densities 
p± = d x if±/2n. Then the effective Hamiltonian reads 

H = r dx [{d^+f + («9^-) 2 ] + Ai cos( v /2^^+(0)), (34) 

where the radiative boundary conditions are 3p + (L) — p + (—L) = 2g + I®/ev + , and 
3p_(L) + p_(— L) = 0. Therefore, only the ip + boson contributes to transport. Next 



we introduce effective fermion fields, ip(x) = y k F /2n exp(i(p + (x)). It is convenient to 
introduce yet another fermion ip(x) via ip(x) = (c + c')^(x), where c is an auxiliary 
fermion. This standard procedure is explained in detail in the textbook [2]. In the 
fermionized framework, Eq. (J34|) becomes 



H = -'^± dx i^cW> + t/^(c + c f ) (HO) - 4'HO)) , (35) 



9+ J V 2 9 

with = irAl/D determining the degree of asymmetry of the junction. In this form, 
the Hamiltonian can be diagonalized immediately, and following Refs. [2 Hj , a self- 
consistency relation determining 51^ follows 

r j / a / \t , (I , AB + 27ri(J 1 °-^ Al /2)/e ^ 

W Al = (eA B /7r)W> ^- + j , (36) 

where ip is the digamma function. The currents Jj are then given by Eq. (J3T|) . with 
obtained as the self-consistent solution of Eq. (j3*6^) . One checks easily that the 
T = corrections for small eAs/If predicted by this solution coincide with our 
previous perturbative results, see Eq. (|3^j) for g + = 1/2. In the opposite limit of strong 
asymmetry, we find the expansion 

= 1 - (7r 2 /6)(A7 e A B ) 2 + . . . (37) 

This corresponds to the highly asymmetric case with almost no current through tube 1 
and the currents through tubes 2 and 3 given by / 2 / 3 = Iy 3 + 1^/2. 

Remarkably, this procedure allows to nonperturbatively treat the full crossover 
between a symmetric Y junction and a single SWNT that is approached by the tip of 
another SWNT at some point in its bulk. The corresponding asymptotic power law 
exponents therefore correspond to the sum of end and bulk exponents in the tunneling 
density of states j2] • Finally, we note that the situation of a perfect beam splitter cannot 
be described by our model. This is due to the initial symmetry of our setup without 
impurities, while a beam splitter can only be realized if transport between two of the 
three tubes is totally suppressed 



5. Conclusions 



In summary, we have proposed a generalization of Landauer-type transport theory to 
iV-terminal starlike interacting nanotubes. Due to the formulation of three different 
boundary conditions, reflecting the coupling to external leads, charge conservation, and 
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density matching conditions at the node, it has been possible to find an analytical 
solution for the nonlinear conductance matrix of the symmetric system for arbitrary 
electron-electron interaction. Based on the symmetric solution, we have constructed 
a wide class of asymmetric transmission matrices by adding impurities close to the 
node in the individual nanotubes. The asymmetric problem can be also be solved 
in a nonperturbative fashion if just one impurity is put into a nanotube arm. For 
simplicity, we have discussed this solution for the special interaction strength g = 3/8, 
corresponding to the Toulouse limit in a folded picture. 

A largely open issue is the question of noise and particle statistics in such a nonchiral 
multi-terminal Luttinger liquid network. For low energies, the disconnected-node fixed 
point seems to hinder the observation of fractional quasiparticle cross-correlations, 
because tunneling of electrons from one tube to the other is the first allowed process 
once the applied voltage is raised. However, for intermediate energies, which could be 
temperature and/or applied voltage in experiments, the disconnected fixed point opens 
up and quasiparticle tunneling between different tubes might give distinct features to 
the current cross correlations and therefore allow to study quasiparticle statistics. These 
questions will be studied in future work. 

Acknowledgments 

We wish to thank H. Grabert, A. Komnik, I. Safi, and C. Texier for useful discussions. 
Support by the DFG under the Gerhard-Hess program and under Grant No. GR 638/19 
is acknowledged. 

References 

[1] S. Datta, Electronic transport in Mesoscopic Systems (Cambridge University Press, 1995). 

[2] A.O. Gogolin, A. A. Nersesyan, and A.M. Tsvclik, Bosonization and Strongly Correlated Systems 

(Cambridge University Press, 1998). 
[3] C. Dekker, Physics Today 52(5) (1999) 22. 

[4] R. Egger and A.O. Gogolin, Phys. Rev. Lett. 79 (1997) 5082; C. Kane, L. Balents, and M.P.A. 

Fisher, Phys. Rev. Lett. 79 (1997) 5086. 
[5] M. Bockrath, D.H. Cobdcn, J.Lu, A.G. Rinzlcr, R.E. Smalley, L. Balents, and P.L. McEuen, 

Nature 397 (1999) 598; Z. Yao, H. Postma, L. Balents, and C. Dekker, Nature 402 (1999) 273; 

H. Postma, T. Tccpcn, Z. Yao, M. Grifoni, and C. Dekker, Science 293 (2001) 76. 
[6] R. Egger and H. Grabert, Phys. Rev. Lett. 77 (1996) 538; 80 (1998) 2255(E); Phys. Rev. B 58 

(1998) 10761. 

[7] R. Egger, H. Grabert, A. Koutouza, H. Saleur, and F. Siano, Phys. Rev. Lett. 84 (2000) 3682; A. 
Koutouza, H. Saleur, and F. Siano, J. Phys. A: Math. Gen. 34 (2001) 5497. 

[8] I. Safi and H.J. Schulz, Phys. Rev. B 52 (1995) 17040; I. Safi, Eur. Phys. J. B 12 (1999) 451. 

[9] D. Maslov and M. Stone, Phys. Rev. B 52 (1995) 5539; D. Maslov, Phys. Rev. B 52 (1995) 14368. 
[10] V.V. Ponomarcnko, Phys. Rev. B 52 (1995) 8666. 
[11] A. Furusaki and N. Nagaosa, Phys. Rev. B 54 (1996) 5239. 

[12] Ya.M. Blantcr, F.W.J. Hckking, and M. Biittiker, Phys. Rev. Lett. 81 (1998) 1925. 
[13] Y. Oreg and A.M. Finkelstcin, Phys. Rev B 54 (1996) 14265. 

[14] A. Kawabata, J. Phys. Soc. Jpn. 65 (1996) 30; A. Shimizu, J. Phys. Soc. Jpn. 65 (1996) 1162. 



Nanotube Y junctions 



15 



[15] A.Yu. Alekseev, V.V. Cheianov, and J. Frohlich, Phys. Rev. B 54 (1996) 17320. 

[16] S. Chen, B. Trauzettel, and R. Egger, Phys. Rev. Lett. 89 (2002) 226404. 

[17] W.Z. Li, J.G. Wen, and Z.F. Ren, Appl. Phys. Lett. 79 (2001) 1879. 

[18] M. Terrones, F. Banhart, N. Grobert, J.C. Charlier, H. Terrones, and P.M. Ajayan, Phys. Rev. 
Lett. 89 (2002) 075505. 

[19] J. Li, C. Papadopoulos and J.M. Xu, Nature 402 (1999) 253; C. Papadopoulos, A. Rakitin, J. Li, 

A.S. Vedeneev, and J.M. Xu, Phys. Rev. Lett. 85 (2000) 3476. 

[20] M.S. Fuhrer et ai., Science 288 (2000) 494; J.W. Janssen, S.G. Lemay, L.P. Kouwenhoven, and C. 

Dekker, Phys. Rev. B 65 (2002) 115423. 

[21] J. Kim et ai, J. Phys. Soc. Jpn. 70 (2001) 1464. 

[22] C. Nayak, M.P.A. Fisher, A.W.W. Ludwig, and H.H. Lin, Phys. Rev. B 59 (1999) 15694. 

[23] I. Safi, P. Devillard, and T. Martin, Phys. Rev. Lett. 86 (2001) 4628. 

[24] A. Crepieux, R. Guyon, P. Devillard, and T. Martin, Phys. Rev. B 67 (2003) 205408. 

[25] S. Lai, S. Rao, and D. Sen, Phys. Rev. B 66 (2002 ) 165327. 

[26] C. Chamon, M. Oshikawa, and I. Affleck, preprint cond-mat/0305121 

[27] J. L. D'Amato, H. M. Pastawski, and J. F. Weisz, Phys. Rev. B 39 (1989) 3554. 

[28] S. Ami and C. Joachim, Phys. Rev. B 65 (2002) 155419. 

[29] A.N. Andriotis, M. Menon, D. Srivastava, and L. Chernozatonskii, Phys. Rev. B 65 (2002) 165416. 

[30] C. Tcxicr and G. Montambaux, J. Phys. A: Math. Gen. 34 (2001) 10307. 

[31] B. Trauzettel, PhD Thesis (University of Freiburg, 2003). 



